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Aims. Our goal is to set constraints on the evolutionary state of the dense core Cha-MMS 1 in the Chamaeleon I molecular cloud. 
Methods. We analyze molecular line observations carried out with the new submillimeter telescope APEX. We look for outflow signatures 
around the dense core and probe its chemical structure, which we compare to predictions of models of gas-phase chemistry. We also use the 
public database of the Spitzer Space Telescope (SST) to compare Cha-MMSl with the two Class protostars IRAM 04191 and L1521F, which 
are at the same distance. 

Results. We measure a large deuterium fractionation for N2H + (11 + 3%), intermediate between the prestellar core L1544 and the very young 
Class protostar LI 52 IF. It is larger than for HCO + (2.5 ± 0.9%), which is probably the result of depletion removing HCO + from the high- 
density inner region. Our CO(3-2) map reveals the presence of a bipolar outflow driven by the Class I protostar Ced 1 10 IRS 4 but we do not 
find evidence for an outflow powered by Cha-MMSl. We also report the detection of Cha-MMSl at 24, 70 and 160 jim by the instrument MIPS 
of the SST, at a level nearly an order of magnitude lower than IRAM 04191 and L1521F. 

Conclusions. Cha-MMSl appears to have already formed a compact object, either the first hydrostatic core at the very end of the prestellar 
phase, or an extremely young protostar that has not yet powered any outflow, at the very beginning of the Class accretion phase. 
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1. Introduction 

The study of the earliest phases of star formation, the so- 
called prestellar and Class stages, is essential to under- 
stand the origin of the st ellar initial mass function (e.g. 
IWard-Thompson et all2 006). As one of the closest active sites 
of low-mass star formation, the Chamaeleon I dark cloud is 
an excellent target to make progress in this field. However, 
its population of prestellar cores is not known and only 1 or 
2 Cl ass protostellar candidates have b een reported so far 
(see Reipurth et al. 1996; Froebrich 2005). The most promis- 
ing of them, Cha-M MSl, is embedded in a C 18 O(l-0) clump 
jHaikala et alJl20o"^> . in the Cederblad 110 region where sev- 
eral young stellar objects have been id entified in the infrared 
(e.g. lPrusti et alJl99lUPersi et all200ll) . Using a tentative far- 
infrared detection, Lehtinen et al. (2001) derived a low temper- 
ature (Tboi = 20 K) and luminosity (Lboi = 0.45 L Q ) for Cha- 
MMS 1 . T hey proposed that it is a Class protostar, in agree- 
ment with iReipurth et aD {l996) who suggested that it could 
be the driving source of the H erbig-Haro obj e cts an d the CO 
outflow seen nearby. However iLehtinen et alJ J2003) failed to 
detect cm-wave emission with the ATCA interferometer, sug- 
gesting it is still prestellar. Here, in an effort to clarify its evolu- 



tionary state, we report molecular line observations carried out 
with the new submillimeter telescope APEX 1 and the infrared 
detection of Cha-MMSl at 24, 70 and 160 fjm by the Spitzer 
Space Telescope. 

2. Observations 

We observed the den se core Cha-MMSl (a 2 ooo = 1 l h 06 m 31 s .7, 
62000 = -77°23'33", iReipurthet alii 19961) in July, September 
and November 2005 with the APEX telescope (Giisten et 
al., this volume). The double-side-band heterodyne receiver 
APEX-2A (Risacher et al., this volume) was tuned to the 
molecular transitions listed in Table [2] The N2H + (4-3) line 
was observed simultaneously with the H2D + line in the 1 GHz 
bandpass. The half-power beamwidth was 18" at 345 GHz, 
and the forward and beam efficiencies used to convert antenna 
temperatures T* into main-beam temperatures T m b were 0.97 
and 0.73, respectively. The double-side-band system temper- 
atures (continuum calibration) ranged from 100 to 360 K in 
T* scale. From spectra taken on different days, we estimate 
the calibration uncertainty better than ~ 15%. The backend 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Fig. 1. a) CO(3-2) integrated intensity observed with APEX 
toward Cha-MMSl (cross). The emission integrated over 
[3.6,5.8] km s is shown in greyscale. The blueshifted, less 
blueshifted and redshifted emissions integrated over [-0.2,1.0], 
[1.2,2.4] and [6.6,9.8] km s _1 are plotted as long-dashed, short- 
dashed and thick contours, respectively. The steps are 3<x, 6<x 
and 3<x with rms cr = 0.13, 0.13 and 0.24 K km s , respec- 
tively. The HPBW is shown on the right. T he star and square 
symbols mark the position of near-infra red JPersi et al.ll200lh 
and 3.5cm sources jLehtinen et alJ2003l) . b) CO(3-2) position- 
velocity diagram along the axis going through IRS 4 and Cha- 
MMS 1 . The dotted line indicates the systemic velocity of Cha- 
MMS 1 . The other lines show the limits of the velocity ranges 
used for the contour maps. The logarithmic temperature scale 
varies from 0.3 to 6 K. 

was a Fast-Fourier- Transform spectrometer with 16384 chan- 
nels and a channel spacing of 61 kHz yielding an effective 
spectral resolution of about 120 kHz (Klein et al., this vol- 
ume). The telescope pointing was checked on R-Dor, 07454- 
7112, IRAS 15194-51 15, or Mars and found to be accurate to 
~ 4" (rms). The telescope focus was optimized on Saturn or 
Mars. The observations were performed in position-switching 
mode using the APECS software (Muders et al., this vol- 
ume). The data were reduced with the CLASS software (see 
http://www.iram.fr/IRAMFR/GILDASl. 

3. Analysis: Kinematics and deuteration 

3.1. No evidence for an outflow driven by Cha-MMS1 

We mapped the region around Cha-MMSl in CO(3-2) with 
a spacing of 15". The spectra are deeply self-absorbed. They 
are very broad with stronger wing emission in the northeast- 
ern part. Figure shows maps of the emission integrated 
over the core and the wings of the line. The blueshifted 
and redshifted emissions are clearly associated with the Class 
I protostar Ced 110 IRS 4 and very likely trace a bipolar 
outflow. Their morphology and the position- velocity diagram 
shown in Fig.[TJ> fall into case 1 of Cab rit & Bertoutl |l990), 
which suggests that the outflow axis is close to the line of 
sight (i S= 30°). The outflow geometry is consistent with 
the geometry of the bipolar nebula seen in the near infrared 
around IRS 4 dZinnecker et alJll999h . The emission integrated 
over less blueshifted velocities traces the same lobe close to 
IRS 4 but in addition shows a prominent extension toward the 



CO 

g CV? 


4 i i i i | i i i i | i 
: HC0 + (4-3) 


1 1 1 1 1 1 1 1 1 1 L 




_ 1 1 | 1 1 1 | 

^N 2 H + (3-2) A 


I I I | I I L 


lO 




uJ^ 

1 , , , , l , , , , l , 


, , 1 , , , , 1 , , , ,- 




: , ' ; '"i " , , , V 


, ' ''! , r i' 1 '' , r 


0.5 


CO 

o 

2" d 


j_i i i i I i i i i I i 
; H i:i C0 + (4-3)y 


1 1 1 1 1 1 1 1 1 1 1 l_ 




i i i i i i i 

7 N 2 H + (4-3) 


I I I I I I L 


m 


CV? 

S ° 
^ o 


~i i i i i i i i i i i i 


I I 1 I I I I 1 I I I ,- 




: , , i T , , i 


i , , ,: 


d 
o 


ID 

g ^ 

J o 
o 


4 i i i i | i i i i | i 
LDC0 + (4-3) 

■? t i i t. I i i i t. I i 


1 1 1 1 1 1 1 1 1 1 
i t I t i i i I t i i r 




i i I i i i I 
7N 2 D + (4-3) 

- , , i , , , i 


! ! ! | ! ! I_ 

r 


0.20.40.6 


o 


_ i i i i I i i i i I i 
- l)CO + (5-4) 


i i I i i i i I i i i i_ 




i i I i i i I 
_ H 2 l) + (l 10 -l n ) 


i i i I i i i_ 


LO 

d 


















" , , , , 1 , , , , 1 , 


: : 1 I I I I 1 I I I 




, , I , , U , I 




o 



23456 2 4 6 

Velocity (km s ) Velocity (km s ) 

Fig. 2. Spectra observed with APEX toward Cha-MMSl (his- 
togram). Synthetic spectra corresponding to the model of 
Sect. l3.3l are superimposed (thick line). The hyperfine structure 
of the N2H + and N2D + lines is shown under each spectrum with 
statistical weighting. The dotted line indicates the systemic ve- 
locity assumed for the modeling. 



near-i nfrared s ource NIR89, a potential Class I young brown 
dwarf tPersi et al.l l2001). This additional wing emission could 
be associated with another outflow driven by NIR89, since 
IPrusti et all dl 99 ll) detected an elongated redshifted emission 
in CO(1-0) on the other side of NIR89 (see their Fig. 5). This 
potential second outflow would have an intermediate inclina- 
tion (20° ^ i 5 70°) since its two lobes do not overlap. On 
the other hand, we do not find significant wing emission in the 
vicinity of Cha-MMSl, which indicates that it does not drive 
an outflow on the scale probed with APEX. 

3.2. Physical structure of Cha-MMS1 

We derive the density structure of Cha-MMSl from 1. 3mm 
conti nuum measurements done with SEST jReipurth et alJ 
Il996l) . Assuming a uniform dust temperature of 10 K, opti- 
cally thin emission and a dust opacity /ci.3 mm = 0.01 cm 2 g _I 
(Ivan derTaketalJ fl999). we derive a mass of 0.40 M in the 
22" SEST be am and 1 .0 M n in the region within the 5<x con- 
tour level of iReipurth et alJ Jl996l) . which has a deconvolved 
diameter of ~ 32". If the dense core has a flat inner region and 
a density decrea sing as r~ 2 , which is t ypical for prestellar con- 
densations (e.g. lBacmann et al.l l2000). then the two measure- 
ments yield a central density bh,, c = 3.3^ ^ x 10 6 crrT 3 and a 
radius of t he flat inner region r^ , - 1800^° AU at a distance 
of 150 pc dKnude & H0g|l998l) . 

3.3. Molecular abundances and deuteration 

The spectra obtained with APEX toward the continuum peak 
position of Cha-MMSl are presented in Fig. |2 and the results 
of single Gaussian and 7-component HFS fits done with the 
CLASS software are written in Tabled The systemic veloci- 
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Table 1. Results of single Gaussian (top) and 7-component 
HFS (bottom) fits to the spectra observed toward the center of 
Cha-MMSl. 
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H 13 CO + (4-3) 
DCO + (4-3) 
DCO + (5-4) 
H 2 D + (l 10 -l,i) 


4.42(8) 
4.42(1) 
4.47(4) 
4.57(4) 


0.53(3) 
0.51(3) 
0.46(8) 
0.63(7) 
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N 2 H + (3-2) 
N 2 H + (4-3) 
N 2 D + (4-3) 


4.44(1) 
4.65(2) 
4.44(2) 


0.42(1) 
0.50(5) 
0.41(3) 


24(2) 
2.3(7) 
0.53(3) 


24(2) 
3.4(14) 
0.5(0) 


50 
81 
44 


0.40(1) 
0.49(5) 
0.40(3) 



The statistical uncertainties in parenthesis are given in units of the 
last digit. means the parameter was fixed. For Vi sl , they include 
the contribution of the frequency uncertainties listed in Table|2| 
The non- thermal linewidth was computed assuming T K = 10 K 
and, for the Gaussian fits, optically thin emission. 



ties derived from the fits agree quite well except for N2H + (4- 
3) and H2D + (lio-ln). Since these two lines were observed si- 
multaneously in the same bandpass, we suspect an instrumen- 
tal problem and have shifted these two spectra by -0.20 MHz 
in Fig. |2 to compensate for this difference. In the following, 
we assume a systemic velocity of 4.43 km s _1 . The opacity 
of N2D + (4-3) is small and hence not well constrained by the 
HFS fit. We fixed it to 0.5 (see model below). The opacity of 
N2H + (3-2) is also uncertain, although this is not seen in the sta- 
tistical uncertainty 2 . Assuming a kinetic temperature of 10 K, 
the non-thermal linewidths are between ~0.4 and ~0.5 km s . 
The linewidths derived from the HFS fits should not be affected 
by optical depth effects. They agree with each other within 2<x 
and we find a weighted mean AVnt = 0.40 + 0.02 km s _1 . This 
corresponds to a non-thermal dispersion cr NX = 0.17 + 0.01 
km s _1 which is on the same order as the thermal velocity dis- 
persion at 10 K. There is therefore a significant amount of non- 
thermal motions in the dense core. 

The HCO + (4-3) spectrum is strongly self-absorbed and 
asymmetric, the red peak being stronger than the blue one, and 
the low-optical depth H 13 CO + (4-3) line peaks in between. The 
asymmetry of the HCO + (4-3) line actually changes across our 
small map: the blue peak is stronger than the red one in the 
southern part belo w (0,-10"). Th is pattern could result from 
rotation (e.g. iRedman et alJl2004h or a combination of infall 
and rotation, provided the actual center of the core is slightly 
offset from our "central" position toward the South (by ~ 5"). 
However we did not detect a clear overall velocity gradient in 
our N2H + (4-3) map over 50" that would unambiguously reveal 
the presence of rotation. 

We used the Mo nte-Carlo -based radiative-transfer code 
MAPYSO (see iBelloche etaD 120021 and references therein) 

2 The self-absorption of the central group of hyperfine components 
lowers its "mean peak temperature" and the opacity derived by the 
HFS fit is then artificially increased. For test purposes, we masked out 
the self-absorption and the opacity was reduced to 19.6 ±1.1. 



Table 2. Gas temperature, non-thermal broadening and abun- 
dance of each molecule for the best-fit model. The peak opacity 
of each modeled transition is also given, as well as an estima- 
tion of the total opacity of the hyperfine multiplets if they had 
been modeled. 
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0.35 
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8.8 


0.35 
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0.50 


2.2(-10) 


15 
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372672.560(13) 
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3.6 


N 2 D + 


4-3 


308422.322(10) 
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0.49 


0.51 


o-H 2 D + 


1 1 o- 1 1 1 


372421.385(10) 
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0.50 
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0.19 





(1) from the CDMS catalog as of Feb. 2006 fsee lMiiller etal . 2005). 
The uncertainties in parenthesis are given in units of the last digit. 

(2) X(p) means Xx 10''. 



to model the emission in HCO + , H 13 CO + , DCO + , N 2 H + and 
N?D + . For all spec ies, we used the collision rates of HCO + with 
H 2 JFlowerll 19991) . We modeled the N 2 H + and N 2 D + lines as 
simple lines, without the hyperfine structure. The H2D + emis- 
sion was modeled with the Monte-Ca rlo program RATRAN 
(see Hogerheiide & van der Takll2000l) . We used the density 
profile described in Sect. 13.21 and the velocity was set to zero. 
We assumed uniform gas temperature, abundance and non- 
thermal broadening for each molecule. We fixed the isotopic 
ratio ^ to 70. 

The synthetic spectra of the best-fit model are overlaid in 
Fig.|2]and the parameters are shown in Table|2] The total opac- 
ities that we estimate from T pea k using the statistical weights 
agree well with the opacities derived from the 7-component 
HFS fits (see Tableland Footnote[5). The non-thermal broad- 
ening derived for HCO + and its isotopologues is more reliable 
than the one measured in Table[T]because the opacity broaden- 
ing is properly taken into account in the modeling. On the other 
hand, it is less reliable for N2H + and N2D + since part of it mim- 
ics the broadening by the hyperfine structure not implemented 
in our model. The non-thermal broadening derived for H2D + is 
larger than for HCO + and its isotopologues, but given the low 
signal-to-noise ratio of the H2D + spectrum, it may not be sig- 
nificant. The one measured for N2H + and N2D + in Table^lies 
in between. From the abundances obtained in Table[2]we derive 
the following deuterium fractionations: [j^+j = 11 ± 3% and 

[hco + ] = ^ ± 0-9%- The deuterium fractionation is therefore 
~ 4 times larger for N 2 H + than for HCO + . 

4. Implications: Evolutionary state of Cha-MMS1 

The deuterium fractionations of HCO + and N2H + derived in 
Sect. I3.3I are 3-4 orders of magnitude larger than the cos- 
mic ratio 7§J ~ 10~ 5 . This enhanced molecular deuteration 

|HJ I II I 

is typical for low-mass dense cores (e.g. [Williams etal. 1998; 
ICraosi et alJ l2005: Paris e et aDl200rj|) and correlates wel l with 
CO depletion in prestellar cores JBacmann et aDl2003l) . It is 
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understood as a result of CO freeze out onto the grain sur- 
face, CO being the major destroyer of H?D + , the key ion in 
the molecular deuterium chemistry ("e.g. ICaselli et alj E0Q3). 
To understand the different deuteration degrees measured for 
HCO + and N 2 H + in Cha-MMSl, we compare our results with 
the predictions of chemical models. Roueff et"al] (|2005) cal- 
culated the steady state of models of gas-phase chemistry at 
10 K. Their three models, Model 1, 2, and 3, have densities 
of 10 4 , 10 5 , and 10 6 cm 3 , and C and O depletion factors of 
1, 5, and 15, respec tively. At 10 K, the ortho-to-para ratio of 
H2D + is close to 1 fcierrich et alJl2002l) . so we deduce a total 
H2D + abundance of ~ 9 x 10", in agreement with Model 3 
and clo se to the value d erived for the evolved prestellar core 
LI 544 ( I van der TalfeOOrJ) . For N 2 H + , both the abundance and 
the deuterium fractionation in Table |3 agree well with the pre- 
dictions of Model 3, while they are close to Model 2 for HCO + . 
The degree of deuterium fractionation in HCO + traces there- 
fore lower densities than the N 2 H + one. This can be under- 
stood if HCO + is depleted at high densities, which is expected 
in dense cores since HCO + is chemically related to CO. To 
test this scenario, we computed a model for HCO + and iso- 
topologues with a hole of radius 1800 AU and abundances 2.5 
times larger: the fit to the observed spectra was as good as for 
the model shown in Fig. [21 suggestin g that HCO + m ay indeed 
be depleted in the inner region. ICrapsi et alJ (12005) derived a 
N 2 H + deuterium fractionation of 1 6-23 % for L 1 544 and 5-10% 
for L1521F, which is now known as a Class protostar since 
its detection by the S pi t zer Space Telescope (SST) in the near 
infrared (e.g. frerebev et alJl2005h . Our value for Cha-MMSl 
lies in between, which suggests an evolutionary state between 
L1544 and L1521F if the |£4£t ratio can be used as a chemical 
clock. 

The higher spatial resolution provided by APEX and 
the higher densities traced by CO(3-2) allowed us to re- 
solv e the blueshifted wing emission see n earlier in CO(l- 
0) JMattilaet alJll989t IPrustiet al1ll99ll) into two outflows 
- the second one being tentative -, and assign their driving 
sour ces. The main out flow is not powered by Cha-MMSl, 
as iReipurth et alJ \ 19961) suggested, but by the Class I pro- 
tostar Ced 110 IRS 4. This is consistent with the non- 
dete ction of Cha-MM Sl in centimeter continuum emission 
jLehtinen et aljE003l) . Since we did not find evidence for an 
outflow driven by Cha-MMSl itself, the dense core must be 
in an earlier evolutionary state than the typical young Class 
pro tostar IRAM 0419 1 which does have a powerful out- 
flow jAndre et alJ Il999l) . However, using the SST archive 
(http://ssc.spitzer.caltech.edu/), we discovered a source de- 
tected at 24, 70 and 160 fim by the instrument MIPS (PID:37, 
REQID:39621 12), offset from Cha-MMSl by 5-7" toward the 
East. Given the large beams of SEST (22") and SST (6-18") 
and their pointing rms accuracies (3" and 1.4", respectively), 
we find it conceivable that this mid-infrared source is em- 
bedded in the Cha-MMSl dense core. We measure fluxes of 
2.5 ± 0.8 and 200 + 100 mJy at 24 and 70 fjm. The 24 p.m 
flux is ~ 7 and 10 times smaller than the fluxes we measured 
for IRAM 04191 and L15 21F using the SST data products of 
the C2D legacy program fevans etalJ Eo03). and the 70 /im 



flux ~ 4 and 3 times smaller 3 . In addition, Cha-MMSl was not 
detected at 8 fim by the instrument IRAC (SST public data, 
PID:37, REQID:3960320), at a 3cr level ~ 9 and 7 times lower 
than the peak fluxes of IRAM 04191 and LI 52 IF The detec- 
tion at 24 fim implies that Cha-MMS 1 alre ady contains a com- 
pact hydrostatic object (e.g. Masunag a & Inutsukall200ol) . but 
the weakness and "redder" color of the 24 and 70 fim fluxes 
and the non-detection at 8 fim sugges t that it is less evolve d 
than IRAM 04191 and L1521F (e.g. lYoung & Evansll2005l) . 
since they are at the same distance. Depending on its inclination 
along the line of sight (see lBoss & Yorkell995HWhitnev et"al] 
2003, for its effects on the mid-infrared fluxes), it could be at 
the stage of the first hydrostatic core, i.e. at the very end of the 
prestellar phase, or have already formed an extremely young 
Class protostar that has not yet powered any outflow. At a 
stage intermediate between L1544 and L1521F, Cha-MMSl 
would then be the first object found so close to the very be- 
ginning of the Class accretion phase. 
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